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Abstract 

This paper describes the calculation of the electromagnetic form factor of the K° meson 
at order p 6 of chiral perturbation theory which is the next-to-leading order correction to the 
well-known p result achieved by Gasser and Leutwyler. On the one hand, at order p 6 the chiral 
expansion contains 1- and 2-loop diagrams which are discussed in detail. Especially a numerical 
procedure for calculating the irreducible 2-loop graphs of the sunset topology is presented. On 
the other hand, the chiral Lagrangian produces a direct coupling of the K° current with 
the electromagnetic field tensor. Due to this coupling one of the unknown parameters of C^ 6 ' 
occurs in the contribution to the K° charge radius. 
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1 Introduction 



Chiral perturbation theory has established itself as a powerful effective theory of low energy strong 
interactions. Based on the symmetry of the underlying QCD, chiral perturbation theory produces 
a systematic low-energy expansion of the observables in this regime. Unfortunately, because of the 
non-renormalizability of the effective theory, higher powers in the energy expansion require higher 
loop Feynman integrals and as input an ever increasing number of renormalization constants. The 
p 4 -Lagrangian involves ten free parameters which were determined in the fundamental papers of 
Gasser and Leutwyler W. For the p 6 Lagrangian there are already more than a hundred of them 
@. Many interesting low energy hadronic amplitudes have now been calculated to order p A and to 
order p 6 |J in the case of SU{2) x SU{2) chiral perturbation theory which involves only one mass 
scale. Recent progress in the calculation of massive two-loop integrals allows now calculations to 
order p 6 in the full SU (3) x SU(3) chiral perturbation theory. With so many unknown parameters, 
however, one may ask about the usefulness of such calculations. Here are a few arguments in favour: 

• In a given class of experiments, such as the electromagnetic and weak form factors of the 
light mesons, only a limited number of renormalization constants enter and relations between 
amplitudes can be tested [|| . 

• The unknown constants enter only polynomially and precision experiments could separate 
the unambiguous predictions (see e.g. the subsequent discussion). 

• Knowledge of the exact low-energy functional form of an amplitude may be important for the 
experimental extraction of low-energy parameters such as charge radii. 

• The results may be used in model calculations which predict the polynomial terms. These 
calculations can then be compared with experiment. 

• The question of convergence of the chiral perturbation theory may be addressed. 

We have embarked on the program of a full p e and two-loop analysis of the various semi-leptonic 
form factors of pions and kaons in full SU(3) x SU (3) chiral perturbation theory. From the point 
of view of practicality it seems prudent to begin such an ambitious program with the simplest 
process possible involving the fewest number of renormalization parameters. Here this would be 
the electromagnetic form factor of the neutral kaon. There exists no tree level contribution to 
C(p 4 ) and therefore the one-loop contributions must be finite. At C(p 6 ), however, a tree level term 
exists which couples the K° directly to the electromagnetic field tensor of the external photon, and 
a new renormalization parameter enters. In this paper we will present the results of our 0(p 6 ) 
calculation of the neutral kaon's electromagnetic form factor. It is defined by the matrix element 
of the electromagnetic current between neutral kaon states: 

{K°,p'\J,\K°,p) = (p + p'),F K °( q 2 ), (1) 

where q = p—p' is the momentum of the photon and stands for the electromagnetic current 
carried by the light quarks which is a linear combination of the chiral currents and V^: 

2 1-1 1 

Jn = 3«7/^- 3 5 7m* = V i + ^ V S- ^ 

We neglect the contributions of the heavy quarks. 
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2 The effective Lagrangian and the diagrams 



Chiral perturbation theory is formulated in terms of an effective Lagrangian involving an increasing 
number of covariant derivatives as well as external sources, 

C eff = £W +£W + £( 6 ) + ... (3) 

where 

£< 2 > = ^ Tr(£» M C/OTt) + L. Tr( X f/t + ^t) (4) 

£ (4) = Li {Tr (D^UD^tf)} 2 + L 2 Tr (Z^C/Z^) Ti{D^UD v U ] ) (5) 

+ l 3 TtiPpUmflDjjirrf) + l 4 TtiD^umj^) Tr( x tf + u x ] ) 

+ L 5 Tr(£> M £/ZWt( x [/t + ^t)) + L(j {Tr( x [/t + f^t)^ 
+ L 7 {Tr( X t[/ - C/t x )} 2 + L 8 Tr( X C/t x[ /t + ^t^t) 
- iL 9 Tr(L,„,mjD , Vi + R^D^D'V) + L w Tr (L^UR^U^) 
£ (6) = -^F^Tr(9^C/t) + ... ; ( 6 ) 

U (x) = exp[i(j)(x) / F] is the unitary 3x3 matrix made up of the Goldstone fields, F is the pion decay 

constant in the chiral limit, the covariant derivative involving external vector- and axial vector 

sources of the corresponding currents, F^ u is the electromagnetic field tensor, and x is related to the 

quark mass matrix. In £( 6 ) we have only displayed the terms relevant to the K° form factor. The 
£(6) 

term is extracted from reference j|| where in their notation f3 = ^F 2 (m? K — m 2 r )(2i?24 — -625)- 
Each term in C e jj produces vertices of 2, 4, 6, . . . mesons with or without attached photon. 
According to Weinberg's power counting theorem || the form factor at order p 6 of the chiral 
expansion is the sum of all 2-loop diagrams with vertices from £( 2 ) , all 1-loop diagrams with vertices 
from £( 2 ) and one vertex from and a tree graph with two vertices from or one vertex from 
£( 6 ). Each dia gram has two external meson lines representing the incoming and outgoing K° and 
one external photon line. All relevant diagrams for the unrenormalized K° form factor are given in 
fig. 1. To indicate the origin of the vertices in the Feynman diagrams we use the following notation: 

£(2) 

-vertices are denoted by filled circles •, £^ 4 ^-vertices by filled squares ■ and an £^ 6 ^-vertex by 
an open square □. Only the vertex is new. 

In each loop diagram, the internal lines symbolize an arbitrary tt, K, or rj meson and they must 
be summed over. Due to the derivative couplings of C e jj the Feynman rules of the vertices are quite 
involved. For the relevant vertices they are discussed in appendix |C| 

The diagrams (la) and (lb) represent the well-known result of order p 4 which is the leading 
order approximation of the K° form factor Q. The pure p e diagrams consist of two groups: 

• the reducible diagrams, i.e. the 1-loop diagrams (2a)-(2e), the 2-loop diagrams (3a)-(3f), 
where the loops are independent of each other, and the tree graph (4) 

• the irreducible 2-loop diagrams (5a)-(5c). 

All diagrams of the first group can be reduced to elementary 1-loop integrals. The results are 
presented in section ||. The calculation of the irreducible 2-loop diagrams is much harder to carry 
out: since in SU(3) chiral perturbation theory there are three different mass scales of the same 
order of magnitude (rn n , mx, Tn„), these 2-loop integrals can no longer be expressed by elementary 
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analytical functions. In section ^| we discuss their contributions to the K° form factor, and in 
appendix [B| we describe our numerical algorithm for the evaluation of these diagrams. 

The contributions of the diagrams can be classified further according to their dependence on 
the momentum q of the photon. Diagrams, where the photon momentum does not flow through, 
are polynomials in q 2 : (lb), (2e), (3d), (5a), and (5c) do not depend on q 2 , (2c) is a polynomial in 
q 2 of degree < 1, and (4) has at most degree 2. Since the normalization of the K° form factor is 
fixed by the charge of the particle due to the Ward-Fradkin-Takahashi identity || , 

F K °(0) = 0, (7) 

the interesting contributions are those with genuine q 2 dependence. Nevertheless, we calculate all 
diagrams, because this procedure offers additional possibilities to check the calculation. 



3 Reducible diagrams 



The diagrams (la)-(3f) are 1-loop integrals or products thereof and can be calculated analytically. 
Inserting the Feynman rules into the vertices yields integrands of the structure 



V(k 2 , kpi , kp 2 ,kl,l 2 , Ipi , lp 2 ) 
P 



(8) 



where p\ and p 2 are the external momenta of the incoming and the outgoing K , k and I are the 
internal loop momenta, the denominator P is a product of scalar propagator factors depending on 
the topology of the diagram, and the numerator V is a polynomial coming from the vertex factors. 
These integrals can be expressed by a set of basic 1-loop functions in the following way: First, p\ 
and P2 are transformed into p = p\ +P2 and the momentum q=p\ — p% of the photon. Then, some 
factors in the numerator are canceled against some propagator factors in the denominator via 



k 



kq 



lq 



P(k, m 2 ) + m 2 



P(k + q, m 2 ) +m 2 - k 2 - q" 



P(l,m 2 ) + m 2 



P(l + q, m 2 ) + m 2 - I 2 - q 2 



(9) 
(10) 

(11) 
(12) 



with P(k, m) = k 2 



m 



What remains are 1-loop 1-point and 2-point functions with tensor numer- 



ators which are decomposable into Lorentz covariants. As our calculation shows, there are in fact 
only two 1-loop functions by means of which all reducible loop contributions of the K° form factor 
can be expressed: On the one hand the scalar 1-point function with its derivative, 



m 



A(m 2 



and A'(m 2 ) 



OA 

dm 2 



(13) 



and on the other hand one specific tensor coefficient of the 2-point function, also with derivative, 

m 

B( q 2 ,m 2 ) = JL -£) and B'(q 2 ,m 2 ) = \£L . (14) 

m 
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These basic 1-loop integrals are denned in detail in appendix [A]. In the following equations we give 
the contribution A of each reducible diagram to the K° form factor in terms of the basic functions 
A and B: 

AM = ±{B(q 2 ,ml)-B(q 2 ,m 2 K )} (15) 
^ = { A ^D ~ A(ml)} (16) 

A (2a) = JL{ B (q 2 ,ml)[2(ml + 2m 2 K )L 4 + 2(ml+m 2 K )L 5 -q 2 L 3 } (17) 

-B(q 2 ,m 2 K )[2(ml + 2m 2 K )U + 4m 2 K L 5 - q 2 L 3 }} 
A( 2fc ) = ^{B(q 2 ,ml)[4(ml + 2m 2 K )L A + 4mlL 5 +q 2 L 9 } (18) 

-B(q 2 ,m 2 K )[A(ml + 2m 2 K )L 4 + 4m 2 K L 5 + q 2 L 9 }} 
A( 2c ) = -j^{A(ml)[4(ml + 2m 2 K )L A +4(ml+m 2 K )L 5 + q 2 L 9 ] (19) 

-A(m 2 K )[4(m 2 n + 2m 2 K )L 4 + 8m 2 K L 5 + q 2 L 9 ]} 

= ™{B(q 2 ,m 2 K )[L 4 (2rn 2 K + ml) + L 5 rn 2 K } (20) 

-£(g 2 , m 2 ) [L 4 (2m 2 + m 2 ) + L 5 m 2 ] 
+£'(g 2 , m 2 K ) m\ [L 4 (2m 2 K + m 2 ) + L 5 m 2 K - 2L % (2m\ + m 2 ) - 2L 8 m|] 
-#V, m 2 ) m 2 n [L 4 (2m 2 K + m 2 ) + L 5 m 2 - 2L 6 (2m| + m 2 ) - 2L 8 m 2 ]} 

A( 2e ) = -l{4 J 4(m 2 )[L 4 (2m| + m 2 )+L 5 m 2 ] (21) 
-4A(m 2 K ) [U(2m 2 K + m 2 ) + L 5 m 2 K ] 

-4A'(m 2 K ) m 2 K [L 4 (2m 2 K + m 2 ) + Lsm 2 . - 2L 6 (2m 2 . + m 2 ) - 2L 8 m|] 
+4A'(ml) ml [L 4 (2m 2 K + ml) + L 5 m 2 - 2L 6 (2m| + m 2 ) - 2L 8 m 2 ]} 

A^ = -^{B(q 2 ,ml) 2 + B(q 2 ,ml)B(q 2 ,m 2 K )-2B(q 2 ,m 2 K ) 2 } 

Am = -^Fl{ 20A (ml)B(q 2 ,ml) + 10A(m 2 K )B(q 2 ,ml) (22) 
+3 A(ml) B(q 2 , m 2 K ) - S0A(m 2 K ) B(q 2 ,m 2 K ) 
-3A(m 2 v )B(q 2 ,m 2 K )} 

A(3C) = ^I^ B ('i 2 ,m 2 K )[A(m 2 v ) + 2A(m 2 K )+A(ml)] (23) 
-4B(q 2 ,ml)[A(ml) + 2A(ml)} 
+2 5'( g 2 , ml) ml[A(m 2 ) - 3 A(m 2 )] 

+ J BV,m!M(m 2 ) (3m 2 +m 2 )} 

A (M) = {3 ^K) iK) + 6 i(m^) 2 - iK) A(m^) - 8 i(m^ 2 (24) 

+2m 2 A' (ml) A(m 2 ) - 6m 2 A' (ml) A(m 2 w ) 
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+ (3m 2 + ml) A! {m 2 K ) A(m^)} 

A( 3e ) = ^{f>A{ml)B( q \ml)+AA{ml)B{q\ml) (25) 
+A(m 2 ) B(q 2 ,ml) - A(ml) B(q 2 ,m 2 K ) 
-8A(m 2 K ) B(q 2 , m 2 K ) - A(m 2 ) B(q 2 , m 2 K )} 

A(3/) = "12^ {^> A (ml) 2 + 3A(m 2 K ) A(ml) + A(m 2 ri ) A(m 2 w ) (26) 
-8A(m 2 K ) 2 - A(m 2 ) A(m 2 K )} 

A«> = ^. (27) 

The contributions contain three masses m n , mx, and m^, but in the isospin limit only two of them 
are independent, m 2 is a short-hand notation for the Gell-Mann-Okubo relation \m\ - \m\ which 
corresponds to the mass of the rj meson in leading order 0(p 2 ). 

Since at order 0(p 2 ) the K° form factor vanishes, the only diagrams which are subject to 
renormalization are the C(p 4 ) diagrams (la) and (lb). Due to renormalization there are three 
modifications of these diagrams: 

• Wave function renormalization of the outer K° lines requires an additional factor 1 + 5Zk 
where SZk is the well-known 0(p 4 ) result (cf. |l]]) 

SZ K = -^{A(m 2 ) + 2A(m 2 K )+A(ml) + 32L±(2m 2 K + ml)+32L 5 m 2 K }. (28) 

• Renormalization of the pion decay constant F w has to be taken into account: the plain 0(p 4 ) 
result contains a factor 1/F 2 where F can no longer be identified with F n at order p e . Instead, 
one has to consider the C(p 4 ) correction F w = F(l + 5f) with (cf. Q) 

S f = ^{M<) + 2Mrnl) + 8L 4 (2m 2 K + ml)+8L 5 ml} . (29) 

• Analogously, the mass renormalizations are to be included, i.e. the masses m appearing in 
(|l5|ff) are related to the physical masses m p h via m 2 h = m 2 + Y,(m 2 h ) where £ stands for the 
0(p 4 ) self energies (cf. |[]]) 

= ^{-A(m 2 v )ml + 3A{ml)ml + A8Uml{2m 2 K + ml) (30) 
+48L 5 m* + 96L 6 m2 (-2m| - ml) - 96L 8 m 4 } 
-Z K (m 2 K ) = {A(m 2 ){3m 2 +ml)+96L 4 m 2 K {2m 2 K + ml) + 96L 5 m 4 K (31) 

+192L 6 m|(-2m| - m 2 n ) - 192L 8 m 4 .} . 

In summary, the total contribution of diagrams (la) and (lb) to the K° form factor at order p e is 
given by 

f 1 + 6Z K ) ( A( la ) + A( lb )) , , . (32) 
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Renormalization corrections of the other diagrams are of order p 8 and can be neglected. 

The sum of all reducible diagrams including the above renormalization corrections takes the 
following compact form: 

F 4 A M; = ±{-A(m 2 v )[4A(m 2 K )-8B(q\m 2 K )-A(ml) + 2B( q 2 ,ml)] (33) 

+ 16 A(m 2 K ) 2 + 10 A{m 2 K ) A(m 2 n ) - 23 A(m 2 n ) 2 - 80 A(m 2 K ) B(q 2 ,m 2 K ) 
+ 32 A(m 2 K ) B(q 2 ,ml) - 28 A{ml) B{q 2 ,m 2 K ) + 70 A{m\) B{q 2 , m 2 ) 
+ 96 B{q 2 ,m 2 K ) 2 - 48 B{q 2 ,m 2 K ) B{q 2 , ml) - 48 B(q 2 , ml) 2 ] 

+ Aq 2 L 3 [B(q 2 ,m 2 K )-B(q 2 ,ml 

+ (4 ml L 5 + q 2 L 9 ) [A{m 2 K ) - 2 B(q 2 ,m 2 K ) - A[m%) + 2 B{q 2 , m\ ) 

Note that the mass derivatives A' and B' and the parameters L4, Lq, and L$ which occurred 
in some diagrams have vanished. The only parameters which are relevant for the K° form 
factor at order p 6 are L3, L5, and Lg. 



4 Irreducible diagrams 

In the irreducible diagrams (5a)— (5c) the 2-loop integrations are not independent of each other as 
they were in the reducible graphs. That is why genuine 2-loop functions enter the stage which 
cannot be expressed by 1-loop integrals only. 

Inserting the Feynman rules yields integrands with a similar structure as in (||), e.g. for diagram 
(5b) 

V(k 2 , kpi,kp 2 , kl, l 2 ,lpi,lp 2 ) 



P(k+p 1 ,mj) P(k+p2,mf) P(k + 1, m|) P(l, mf) ' 
where V is a polynomial of degree equal to the number of vertices. After canceling factors via 



kpi 


1 r 
2~l J 


kp2 


1 r 


kl 


1 r 
2l J 


I 2 


= P(l 



P(k+pi, m\) + m,0 — k 
P(k + p2, m\) + — k 2 



2 2 
Pi 



vl 



c 2 - I 2 



P(l,m 2 ) + m 2 3 



(34) 

1 

(35) 

(36) 

(37) 
(38) 



we are left with reducible integrals which can be calculated analytically, and with some genuine 
2-loop integrals of the sunset-topology, i.e. the 3-point functions 



p / 2 2 2 2 2 2 2\ 

^a ua2 Al '^Pi^P2;'m ,m l ,m2,m 3 ) 



-2D 



2 2 2 2^ 

m 1 ,m 2 ,m 3 " 
d D k d D l 



{lpi) ai {IP2T 2 (k 2 f 



i{2-n) D i{2ir) D P(k + Pl ,ml) P{k+p2,ml) P(k + l,ml) P(l,ml) 
and the 2-point functions 

dPk d D l {lp) a (k 2 ) 13 



S a ,p{P ; m i, m 2, m 3) 



8-2D 



i(2ir) D i(2ir) D P(k+p,mf) P(k + l,m 2 ) P{l,m 



(39) 



(40) 
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In diagram (5b), five different mass flows of intermediate mesons must be regarded: 

A (56) _ a(56) , a(56) A (56) A (56) A (56) {m 

where A^) means that meson r couples to the photon and the other two lines are mesons of type 
s and t. Each mass flow is handled separately, and its contribution to the K° form factor can be 
expressed in terms of the basic 1- and 2-loop functions A, B, S a ,/3, T ai ^ a2: p, where for the latter at 
most the tensor indices 

S 2 ,o, 5i,i> ^1,0) Sb,0; -Sb,!, Sq, 2 , 7b,o,o> ^b,o,i' ^0,0,2, ^0,0,3, ^i i0 ,o, ^1,0,1, ^1,0,2, Ti,i,o, ^1,1,1 



are needed. In the following table we give the detailed result for each mass flow. We omit the 
arguments of the 2-loop functions S a $ and T QljQ2j p: for mass flow A^ they are understood to be 

S a ,p(m 2 K ;m 2 ,m 2 , m 2 , ) and T aua2 ^(q 2 ;m 2 K ,m 2 K ; m 2 , m 2 ,m 2 ,m 2 ). 



A ( 5b ) - - / (421 

K ™ ~ %F±m 2 (Am 2 -q 2 )\ ( ] 

- Ti.1,1 • 32 m 2 K + T 1;lj0 • 16 m 2 K (4m 2 K - q 2 ) + T lfi>2 ■ 16 m 2 K 

+ Ti i0 ,i • 8 m 2 K (-Am 2 K + q 2 ) - T 0i0 ,3 ■ 2 m 2 K + T 0i0 ,2 ■ m 2 K (4m 2 K - q 2 ) 
+ S 2fi ■ 16 (2m| - q 2 ) + • 8 (-4m| + g 2 ) + S lfi • 8 m|(4m| - q 2 ) 
+ S ,2 ■ (6m| - g 2 ) + 5 0) i ■ 2 m|(-4m| + g 2 ) 
+ A(ml) 2 ■ (Am 4 - 8m|m 2 - m 2 K q 2 + 2m 2 <j 2 ) 
+ A(ml) B(q 2 ,m 2 K ) • 4 m|(4m| - q 2 )} 

A ( 5b ) = - / (43) 

24F 4 m| ( 4m 2 - ,2) \ ^ 

- T ,o,3 • 18 ml + T 0i o,2 ■ 3 m|(28m| - 3g 2 ) + T , ,i • 8 m 4 K (-16m 2 K + 3q 2 ) 

+ T ,o,o • 16 m|(4m 2 . - g 2 ) + S , 2 ■ 3 (10m 2 . - q 2 ) + 5 0) i ■ 2 m 2 K {-Um 2 K + 5q 2 ) 
+ 5 ,o • 8 m 4 K (8m 2 K - q 2 ) 

+ A(ml ) A(m 2 ) • (4m| - %m 2 K m\ - m 2 K q 2 + 2m 2 q 2 ) 
+ A(m 2 ) B(q 2 ,ml) ■ 10 m^m 2 , - q 2 ) 
+ A(m 2 ) £(g 2 , m 2 ) • 10 m|(4m| - g 2 )} 

A ™ K = 12F4m|(4m 2 : -g2 ) { ( 44 ) 

- Ti,i,i • 24 m| + T 1)lj0 • 12 m|(4m| - g 2 ) + T 1)0 ,2 ■ 24 m 2 K 

+ r 1A i • 12 m 2 K {-Am 2 K + g 2 ) - T , ,2 • 24 m| + T , ,i ■ 12 m|(6m| - q 2 ) 

+ T ,o,o • 12 m|(-4m| + g 2 ) + S 2fi ■ 12 (2m| - q 2 ) + • 10 (-4m| + g 2 ) 

+ 5 1)0 • 8 m^m 2 , - g 2 ) + ,So,2 • (4m 2 K - q 2 ) + S 0) i ■ 24 m 4 

+ S 0fi -4m 4 {-Wm 2 K + q 2 ) 

+ A(m|) 2 • 3 mli-Aml + g 2 ) 

+ A(m 2 K ) B(q 2 , m 2 K ) ■ 2 m 2 K (-4m 2 K + g 2 )} 

A ( 5b ) = I / (45) 

WK7T 24F 4 m|(4m2 - q 2 ) I 1 J 
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• 72 m 2 K + Ti.1,0 ■ 36 m|(-2m| - 2m 2 + g 2 ) 
+ Ti.0,2 • 24 m| + Ti )0) i ■ 12 m 2 K (-4ml - 4m 2 + q 2 ) 
+ T 1A0 • 12 m 2 K (2m% + 4m|m 2 - m 2 K q 2 + 2m* - m 2 g 2 ) 

- To i0 ,3 • 6 m 2 . + T 0i0 ,2 ■ 3 m 2 K (6m 2 K + 6m 2 - g 2 ) 

+ ^o,o,i • 6 m 2 K (-?>m% - <om 2 K m\ + m 2 .^ 2 - 3m£ + m 2 g 2 ) 

+ T ,o,o • 3 m 2 K (2m 2 K + 2m 2 - g 2 )(m 2 + m|) 2 

+ ^2,0 • 36 (-2m 2 . + q 2 ) + • 2 (-4m 2 . + g 2 ) 

+ 5i )0 • 2 (4m* + 4m|m 2 - m 2 K q 2 - m 2 g 2 ) + S 0>2 ■ 2 (7m 2 . - q 2 ) 

+ 5 ,i • 4 (-7m£ - 7m 2 .m 2 + m 2 K q 2 + m 2 <? 2 ) 

+ 5 ,0-2(7m 2 f -g 2 )(m 2 +m|) 2 

+ A(m 2 K ) B(q 2 , ml) • 22 • m|(-4m| + q 2 ) 

+ A{ml) B(q 2 ,m 2 n ) ■ 8 m|(4m| - q 2 )} 

A ^ = 24F^m 2 (4m 2 - q 2 ) i (46) 
Ti.1,1 • 72 ml + Ti.1,0 ■ 36 m|(-2m| - 2m 2 + q 2 ) 

- Ti.0,2 • 72 m| + Ti )0) i ■ 12 m|(4m| + 12m 2 - 3g 2 ) 

+ Ti i0 ,o • 12 m 2 K (2m K - 4m 2 .m 2 - m|^ 2 - 6m^ + 3m 2 g 2 ) 

+ T , ,3 • 18 m| + T ,o,2 • 3 m|(-2m| - 18m 2 + 3q 2 ) 

+ T 0i o,i • 2 m 2 K (-hm\ + 6m|m 2 + 3m|g 2 + 27m* - 9m 2 g 2 ) 

+ T ,o,o • m|(g 2 - 2m| - 2m 2 ) (-3m 2 + m|) 2 

+ S 2 , • 36 (-2ml + <Z 2 ) + • 30 (4m| - q 2 ) 

+ Si fl ■ 6 (-4m* - 20m|m 2 + m|g 2 + 5m 2 g 2 ) + S 0>2 ■ 6 (-7m 2 . + q 2 ) 

+ 5 0) i ■ 4 (m 4 . + 21m|m 2 - m|g 2 - 3m 2 g 2 ) 

+ So, -2 (m 2 K - 3m 2 )(3m| + 7m 2 .m 2 - m|g 2 - m 2 g 2 ) 

+ A(m 2 K ) A(m 2 ) ■ 4 (12m£ - 8m|m 2 - 3m|g 2 + 2m 2 g 2 ) 

+ A(ml) B(q 2 , m 2 ) • 2 m 2 f (4m| - <? 2 ) 

+ A(m 2 v ) B(q 2 ,ml) • 16 m 2 K (-4m 2 K + q 2 )} . (47) 

Diagrams (5a) and (5c) yield the same contribution because of time reversal invariance. In these 
diagrams, there are three mass flows 

A<» = Affi + A£> + A0*> , (48) 

which can be expressed by the basic functions A and S a ^: 

= 24F±m 2 { ~ 8 + SlA + 12 m ' Sl >° + 5 °' 2 ~ 6 m ' 5011 (49) 
+ 3 m| A(m|) ^(m 2 ) + 2 (2m| - m 2 ) A(m 2 ) 2 } 

= 40 f4 2 { 6 ~ 8 m l ■S'l.o + 3 5*0,2 + 2 m 2 5o,i - 8 m\ 5*0,0 (50) 
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- 6 m% A(m 2 K ) A(m 2 r ) - 4 m| A(m 2 K ) A(m%) 
-2(m 2 K + ml)A(m 2 7r )A(m 2 ri )} 

^kkI = 24F 4 m2 { ~ 4 + 2 Si.i - S , 2 + 3 m 2 5 ,i + m 2 , ^(m|) 2 } . (51) 

Except for special kinematic situations the genuine 2-loop integrals S a ,/3 and T ai)Q2j/ g cannot be 
calculated analytically. In appendix [B] we describe the method how we calculated them by splitting 
them into one part which contains the divergence and can be evaluated analytically, and a second 
part which is finite and can be done numerically. 



5 Evaluation of the diagrams and checks 

In the previous sections all 0(p e ) contributions to the K form factor have been written in terms of 
some basic 1- and 2-loop functions. Before evaluating the contibutions explicitly, we must specify 
a renormalization scheme. 

In our calculation we are using dimensional regularization and the so called GL-scheme which 
is defined in the following way: Each diagram of order 0(p 2n ) is multiplied with a factor e ( 1-n ) a ( £ ) 
where D = A~ 2e is the dimension of space-time and a(e) is given by 

(47r) £ r(-l + £ ) = , (52) 



that is 



vr 2 In ^ 



a(e) = e {l- 1 + \ogATT)+e\- + -)+0{e i ). (53) 



Because of a(0) = the total 0(p e ) result is unchanged in D = 4 dimensions. The reason for this 
modification of each diagram is to eliminate the geometric factor (4tt) £ T(— 1 + e) appearing in the 
1-loop integrals A and B. This renormalization scheme is very similar to the well-known MS-bar 
scheme: the only difference is that in MS-bar the left-hand side of the defining equation ( |52|) is 
- {Air) 6 r(e) so that there a geometric factor of (4-7r) £ T(e) is eliminated from the 1-loop integrals. 

The GL-scheme extends the usual 1-loop scheme introduced by Gasser and Leutwyler BJ in a 
natural way which can be understood from an inspection of the renormalization constants Li of 
£( 4 ) : In D-dimensional space-time they have dimension D — 4 and their dimension is generated by 
the mass scale \i of dimensional regularisation: 

U = ijP-^L^D). (54) 

Li(/j,,D) has the same /i-dependence as a 1-loop integral, because Li itself is independent of fi. It 
can be expanded in a Laurent series around e = in the same way as a 1-loop integral: 



r(-l) 

LfaD) = ^— + Lf\^+eL^^) + 0{e 2 ). (55) 



In the usual 1-loop scheme one chooses 



L) ' = (56) 
LfV) = Lr(/x)-^[l-7 + log(47r)], (57) 
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where are numbers which can be found in O . The second term in L^ ' is constructed so that it 
cancels in the £°-coefficient after multiplication with e~ a ^: 

r(-i) 

L? L (fi,D) := e-^LifaD) = -A— + L f^ GL {^) + e L^ GL {^) + 0{e 2 ) (58) 

with L^ GL ( M ) = Lf"^). 

The dimension of the £( 6 )-parameter /3 appearing in (|6|) can be treated in the same way: 

P = fi 2D - 8 P(fi,D), (59) 

where /?(//, -D) behaves like a 2-loop integral. Its Laurent series in the above GL-scheme is given by 



a(-2),GL at-V),GL(..\ 



e 2 



(60) 



Before discussing the numerical results we enumerate the checks which we performed on our 
calculation: 

• Because of charge conservation, the on-shell current matrix element in (|l|) must not contain 
a term proportional to (p — p')u- This is manifestly the case for each individual diagram after 
inserting the Feynman rules into the vertices. 

• According to the Ward-Fradkin-Takahashi identity B the form factor must be equal to the 
charge of the particle for zero momentum transfer. This is not the case for each individual 
diagram, but for the sum of all reducible and the sum of all irreducible diagrams separately. 

• The most important check of our calculation follows from an analysis of the divergent parts 
of the loop diagrams: their sum must be equal to the negative divergent part of the tree 
graph (4), so that the sum of all p 6 diagrams is finite. Since diagram (4) is a tree graph, it 
is a polynomial in masses and momenta and cannot produce logarithms thereof. Thus, all 
logarithmic terms must cancel in the sum of the divergent parts of the loop diagrams. For 
the irreducible diagrams, we find in the GL-scheme 

« = ^iSfet 1 + 432(4W7{ " 16 -J< 10 ^ " ^(-i) (61) 

-ml (32m| - 80m 2 - 9q 2 ) Lg(m 2 n ) - 3g 2 (4m| - ml) Lg(m 2 ) 
-b {q 2 ;m 2 K , m 2 K ) (160m* - 112m|m 2 - 76m 2 K q 2 + 28m 2 < ? 2 + 9g 4 ) 
-b (q 2 ;ml, m 2 n ) (32m 2 K m 2 w - 8m 2 K q 2 - 80m 4 + 56m^ 2 - 9q 4 ) 
+80g 2 (m|-m 2 )} 

using Lg(m 2 ) = log(m 2 /(47r^ 2 )) +7, and for the reducible 1- and 2-loop diagrams 



red. loop \ irr. loop V , 2 _ 2 \ t ren 

(4vr) 2 F 4 e 



.red. loop A irr. loop ,2 2 wren 

^div.part — div. part ' ( a„\1t?A^\^k 7rj 3 • \ 0Z ) 



In fact, the logarithmic terms vanish in the sum of all diagrams, but no longer separately for 
the group of reducible and irreducible diagrams. 
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The previous two formulae (pT|) and (^) , together with equation (|27|) , fix the divergent part of 
the relevant constant (3 which occurs in the K° form factor: 

2 2 

f-^div — -^3 ' l M J 

Note that there is no 1/e 2 contribution, and as a consequence the given 1/e result is independent 
of the renormalization scheme and the renormalization mass /i, because it is the leading order term 
in the Laurent expansion of (3 w.r.t. e. 

For the numerical evaluation we choose the following input data: 

m K = 495 MeV (64) 
= 0.28m K = 138.6 MeV (65) 
F w = 92.4MeV. (66) 

Note that we have used m 2 as a short-hand notation for the Gell-Mann-Okubo term \vr? K - ^rn 2 
in all contributions listed in the previous sections. The loop integrals are calculated at mass scale 
/x = m p = 770 MeV. At that scale the relevant parameters are given by 

Ll en {m p ) = (-35 ±13) • 10" 4 (67) 
L r 5 en (m p ) = (14±5)-10" 4 (68) 
L r 9 en {m p ) = (69±7)-HT 4 . (69) 

Finally, at 0(p 6 ) not only the ^-coefficient L,f > '' GL = Ll en occurs in the result, but also the coefficient 
L^' GL of the subsequent order e 1 of the Laurent expansion (^|). It can be proved that these 
constants L^' GL always enter the result in a specific combination with some 

£(6) 

parameters, so 

that they are no new degrees of freedom and can be chosen arbitrarily: a change in their values 
would only lead to a redefinition of some £^ constants. For our calculation, we choose 

L^' GL (m p ) = 0. (70) 

With these input data and definitions we obtain the contributions shown in fig. 2 for the 
reducible and the irreducible loop diagrams. The corrections are to be compared to the leading 
order result which is also shown in fig. 2. 



6 Discussion and Results 

Since at order p e there exists a direct coupling of the K° to the external photon, the slope of the 
form factor at zero momentum transfer is not predicted. Diverse models make predictions for the 
direct coupling (3 of Eq. (0), but ultimately only a precise measurement of the K° charge radius 
would fix the constant (3. On the other hand, to extract the neutral kaon's charge radius from data 
a knowledge of it's precise functional dependece on the momentum transfer q 2 is essential. Our 
explicite results come here extremely handy, as only one parameter (the charge radius) needs to be 
fitted to the data. 

Alternatively one may consider the observable F K ° (q 2 )/q 2 in the comparison of theory with 
experiment. Here the unknown constant (3 only influences the vertical position of the predicted 
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Figure 2: Contributions to the K° form factor: the solid line is the leading order 1-loop result, the dotted curves 
denote the corrections at order p 6 due to the reducible resp. irreducible loop diagrams in the GL-scheme (on the 
one hand diagrams (2a)-(3f) and the p 6 contributions of the renormalized leading order diagrams (la)-(lb), on the 
other hand diagrams (5a)-(5c)). The contribution of diagram (4) is missing: It would add a straight line through 
the origin with a slope given by the £^ 6 - ) -parameter (3. 

form factor curve, but not its shape. The predictions are plotted in fig. 3. We do not include the 
errors due to the constants as their effect tends to be proportional to q 2 , and therefore cancels 
in the plot of fig. 3. 

It can be observed that the 0(p e ) correction to the 0(p 4 ) approximation is unexpectedly large. 
This result may indicate a breakdown of chiral perturbation theory or, as the 0(p e ) term represents 
only the first correction to the leading 0{p i ) term, it may be due to the fact that the 0(p 4 ) term 
is accidentally small. A detailed measurement of the form factor of the neutal kaon may provide 
unique information on the convergence of SU(3) x SU{3) chiral pertubation theory. This is because 
the prediction of the shape of F K ° ' (q 2 )/q 2 , in contrast to most other 0{p G ) results, does not involve 
any of the unknown 143 freely adjustable parameters of ||. There exist plans for experiments on 
the form factor of the neutral kaon at CEBAF and MAMY C, and we can only hope that these 
experiments are realized. 

A The relevant 1-loop integrals 

All reducible diagrams contributing to the K° form factor were reduced to two basic integrals in 
section |3[ Here, we give the definitions and the results of these basic integrals in dimensional 
regularisation. 

The first basic 1-loop integral is the 1-point function 



where /i is the scale of dimensional regularisation and D = 4 — 2e is the dimension of space-time. 
What is relevant in a 2-loop calculation are the coefficients of the Laurent expansion of the integral 




(71) 
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Figure 3: K° form factor divided by q 2 , in leading and next-to-leading order of ChPT. The curves are shifted by a 
constant so that they go through the origin. This constant is not predicted by ChPT at next-to-leading order, only 

at leading order. 



around e = up to order e 1 . For A we have 



A{m 2 ) = ^-i + a0 + ea 1 + O(e 2 ) (72) 



with 

(47T) 



2 

m 



2 

m 

ao = (4^ 
m' 

ai ■- 



(4vr) 2 { 1 - Lg(m2) } (74) 

m2 (-, 71-2 T / 2^ Lg(m 2 ) 2 i 

{1 + --^) + -!^} (75) 



where we have introduced the abbreviation 

2 

Lg(?7l 2 ) = log ) + lEuler • (76) 

The second basic 1-loop integral, the 2-point function B(q 2 ,m 2 ), is defined as the coefficient of the 
Lorentz-decomposition of the tensor integral 

a^f-M ^ (77) 

M J i(2ir) D [{k + q) 2 — m 2 } [k 2 — m 2 ] { ' 

coming along with the metric tensor g^ v where the other Lorentz-covariant is q^q u . It can be 
expressed further by A and the scalar 2-point function 

B (q 2 ,m 2 ) = J UU [{k - - 2 _ -L [k2 _ m2] (78) 
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in the following way: Let 

B (q 2 ,m 2 ) = h -^ + b + eb 1 + O{e 2 ) (79) 

B(M _ (a=l + (B) 0+£ ( B ) 1+O (a (80) 

be those parts of the Laurent series which are relevant in a 2-loop calculation. 6_i, 60 and b\ are 
functions of to 2 and q 2 and are given below. The relationship between B on the one side and the 
scalar integrals A and Bq on the other side follows from a tensor decomposition: 

B(q 2 ,m 2 ) = -^—^{2A(m 2 ) + (4m 2 -q 2 )B (q 2 ,m 2 )}. (81) 

The relevant Laurent coefficients of B are given in terms of those of A and Bq via the polynomial 

1 / 2 4 n 

P(xi,x 2 ,x 3 ) = — (x 1 + -x 2 + -x 3 J (82) 

in the following way: 

(£)_i = 2P(a_i,0,0) + (4m 2 - g 2 )P(6_i, 0, 0) (83) 

(B)o = 2P(a ,a_ 1 ,0) + (4m 2 - (/ 2 )P(6 ,6_ 1 ,0) (84) 

(B)i = 2P(ai,ao,a_i) + (4m 2 - q 2 )P{h, 60,6-1). (85) 

The basic Laurent coefficients 6_i, 6 and 61 of the scalar 2-point function P involve logarithms 
and dilogarithms: 

(4vr) 2 6_i = 1 (86) 
(4vr) 2 6 = 2 -Lg(m 2 )-rT 1 (87) 
r 2 
12 

where we have defined 



(47r) 2 6! = 2 + ^ (1 - 2r) + i (Lg(m 2 ) - 2) 2 + r(TiLg(m 2 ) - 27\ - ^T 2 + T 2 + 2T 3 )(88) 



r = ^-f (89) 
Pi = lo g(^j) (9°) 

T2 = log2 ( I 2r") (91) 

Ts = Li2 ( I 2T L ) • (92) 

An infinitesimal negative imaginary part of all masses is understood. 

The Taylor expansions of the Laurent coefficients of B w.r.t. q 2 is needed if one is interested 
only in the small q 2 behaviour of the form factor: 

(4^) 2 (P) = -^!{ L g(m 2 )-l} + ^Lg(m 2 ) + 0(g 4 ) (93) 

(4*) 2 (*)i = ^{(Lg(m 2 )-l) 2 + l + ^}-|!{LgK) 2 + ^} + 0(, 4 ). (94) 

With the help of the above formulae the basic 1-loop functions A and B and therefore all reducible 
loop contributions to the K° form factor are reduced to logarithms and dilogarithms. 
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B The irreducible 2-loop integrals 



In section ^| the irreducible 2-loop diagrams of the if °-form factor were expressed by a set of basic 
2-loop integrals of the 2-point and 3-point sunset topologies S a ^ and T ai a2 p, cf. ( 39 Ff ) . Fig. 4 
shows the diagrams and the flows of momenta. 




k+pi,m 



k+p 2 ,mi 



Figure 4 



Z,m 3 




k +p,mi 



In the case of three different masses (m n , m K , and m„) it is no longer possible to reduce them to 
elementary analytical functions. Therefore, we chose a numerical approach for our calculation which 
will be outlined in the following. The procedure is explained in the example of the 3-point function 
T ai a 2 ,/3'i the- 2-point function S a ^ is simpler and can be treated along the same lines. Tensor 
integrals of the relevant topologies are decomposed into Lorentz covariants and the algorithm is 
applied to the coefficient functions. 

The idea is to split the integral T aija2j g into two parts, 



t ai,a2,/3 



f ai,af2,/3 



+ T, 



Ql,Q2,/3 ' 



(95) 



where a2 g contains the divergences, but is of a simpler structure so that it can be calculated 
analytically, and Q2 g is finite in D = 4 dimensions and can be evaluated numerically. 



Our method of achieving a suitable decomposition ( |9q ) is based on the well-known BPHZ 
regularization procedure @ which makes use of the following property of Feynman integrals: each 
Feynman diagram without subdivergences behaves asymptotically as a polynomial in its external 
momenta. Therefore, an UV-divergence of a subdivergence-free diagram may be extracted by 
subtracting a Taylor polynomial of sufficient degree w.r.t. the external momenta. 

We apply this subtraction of Taylor polynomials in two steps: First, we consider the l- 
subdiagram which is a 2-point function with external momentum k, the loop momentum of the 
other loop integration. Subtraction of a Taylor polynimial w.r.t. k renders the /-subgraph finite 
and yields an integral 



-ai,ct2,/3 



dk 



(k 



2\a 



P, 



k+p 1 ,ttIq ,m\ 



(l - Taylor J J 



dl 



(lpi) ai (lP2 



P 



fc+i,m|P?,' 



(96) 



which has no subdivergences. The diagrammatic structure of T ai a2 p is shown in figure 5. 

In the second step we subtract a Taylor polynomial from T aijCe2 p w.r.t. its external momenta p\ 
and p2- since T aii0l2i g has no subdivergences left, we end up with a finite integral which is defined 
to be a2 g in the decomposition (|95"|). 

This decomposition has the desired properties: On the one hand, the part T£ a2 g containing 
the divergences has the diagrammatic structure shown in fig. 6 and can be reduced to well-known 
analytic functions. This is obvious for the first and third topology of fig. 6 which are products of 
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Figure 6 Propagators marked by a x sign are potentially taken to some integer power. 



1-loop integrals. The second topology in fig. 6 has the structure of a scalar 2-loop vacuum bubble 
which are discussed in detail in ||. 

On the other hand, the finite part a2 a is accessible to numeric evaluation, because seven 
of its eight integrations can be done analytically, and the integrand of the final integration is 
a smooth function consisting of logarithms and dilogarithms. This can be understood from the 
following observation: since the /-subintegral is a 2-point function with external momentum k, it 
may be represented by a dispersion integral 

Ar /aia 2 (C»"l2» m 3»Pl»P2) (Qffs 

aQ , (97j 

(m,2+m,3) 2 

where the k dependence occurs only in the dispersion denominator. Interchanging the k and £ 
integrations we end up with an integral 

d( f ai a 2 (C,m 2 ,ml,pi,p 2 ) / dk — — -r- , (98) 



?(C-fc 2 )' 



where the k integration can be done analytically if one recognizes the dispersion denominator 
£ - k 2 as a propagator with mass £ and momentum k. Thus, the k sub integration is a 1-loop 
3-point function which is given in terms of dilogarithms in ||. 

This procedure of reducing the numeric part Q2 ^ to a 1-dimensional integration is illus- 
trated diagramatically in figure 7: writing the /-subgraph as a dispersion integral and subsequently 
interchanging the integrations has the effect of shrinking the /-subgraph to a single propagator of 
mass £ which must finally be integrated over. 
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The above method of calculating the irreducible 2-loop diagrams is described in more detail 
in [|nj. We implemented the algorithm in a set of REDUCE programs which can be supplied on 
request. 



C Feynman rules 

In this section we discuss the momentum space Feynman rules which are needed for the relevant 
vertices occurring in the diagrams of the if°-form factor (cf. fig. 1). There are ten different vertex 
types which are to be considered: vertices from £,( 2 \ and with an even number of meson 
legs and with or without an additional photon (fig. 8). 






(a) (b) (c) (d) (e) 






(h) (i) (j) 

Figure 8: The vertices which are needed for the K° form factor up to order p 6 . A filled circle (•) denotes a vertex 
from £ (2) , a filled square (■) a vertex from and an open square (□) a vertex from . 

In general, the Feynman rules for a vertex with n meson legs is obtained from C e jj by first 
determining all monoms in the expansion of C e jj which contain exactly n factors eft or derivatives 
thereof. Then, transformation into momentum space is done by replacing each derivative by 
—ip^ where p^ is the momentum flowing into the corresponding leg of the vertex. Finally, the 
vertex must be symmetrized in all its meson fields due to Bose symmetry. 

For the vertices (a) and (b) of fig. 8 the relevant monoms of the expansion of £( 2 > are 

4— = 2^2^([^W^) + Tr( X 4 ) (99) 

4 2 ™ = ^^{lWWV4T:^ (100) 

^Tr(y0 6 ) . 

1440F 4 KXV ' 
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In our calculation we avoided writing down the Feynman rules for a general particle flow in terms 
of the SU(3) structure constants f a bc and d aoc . Instead, we calculated the traces in the Feynman 
rules for each concrete flow of particles separately and inserted the Gell-Mann matrices A a which 
occur in 



= 



a=l 



( tt° + ^v v / 2vr+ V2K+ \ 

\ -fa ) 



(101) 



directly into our REDUCE programs where also the symmetrisation is done. In this way, the 
Feynman rules can be read off from the relevant monoms of C e jj. The only step which requires 
some effort is the determination of these monoms. 

The mass matrix can be expressed in terms of the unrenormalized meson masses and takes the 
following form in the isospin limit: 

X = dmg(ml,ml,2m 2 K -ml) . (102) 

For vertices (c) and (d) of fig. 8 the monoms defining the Feynman rules are 

4 4 !™ = ^ Tr(a M dty Tr( X ) + ^ Tr( X ffty - ^? Tr( X 2 ) Tr( X ) (103) 

" fH Tr (^)} 2 " ^J{^(xV) + Tr( X ^)} , 
4— = ^{lV(^^)} 2 + ||Tr(^^)Tr(^5» (104) 

- 7^{M[<P, d u <f>]cf> V<f>) Tr( X ) + 3 Tr(cW> d^) Tr( X 2 )} 
L 5 



Gp4 -{2 Tr( X 2 d v <j> d»<j>) + 3 Tr( X d u cf> d»cf> </>) - Tt( X 4> 4> d'cf) 
+ Tr( X d v cp <f> 2 d v <\>) - Tr( X drf <f> d u <j) <f>) + 2 Tr(x d v <j) 8 u <j) (j?)) 
+ ^{Tr(x<A 4 ) Tr(x) + 3 [Tr( X 2 )] 2 } + ^ Tr( X ^ 3 ) Tr(^) 



+ 



6F 4 



{ Tr(xV) + 2 Tr( X 0X<A 3 ) + 3 T^X^W) + 2 Tr(x</> 3 X</>) } • 



The remaining vertices (e)-(j) contain an even number of mesons and one photon line. In the 
following we give the corresponding Feynman rules for an arbitrary external boson, e.g. a photon 
or a W ± . The external boson enters the Lagrangian C e g as a gauge field of the chiral symmetry 
group in the covariant derivative 

D„U = d^U + iUl^-ir^U (105) 
and can be expressed by the Gell-Mann matrices T a = A a /2: 

l, = Y,llT\ r , = j2r;T a . (106) 

a=l a=l 
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In case of electromagnetic interaction, the gauge boson is the photon and is given by 

h = = -eA^ diag(|,-i,-i) . (107) 
A left-handed gauge boson l a couples to the meson current 



(108) 



the analogue is true for a right handed one. J^/a R is therefore the relevant monomial part of 
C e ff which yields the Feynman rules for the meson vertices with one external boson. Since C e jj is 
symmetric under intrinsic parity transformation, and J^ a are not independent of each other: 

J»A U ) = J U U] ) reS P- J ZaW = J M,a(-^) ■ ( 109 ) 

Therefore, it suffices to specify only jfz a for the relevant vertices. 

The Feynman rule of the vertex from with n mesons and one external boson is given by 

■ rp2 -n n—l 
fc=0 

Vertices (e)-(g) are special cases of this formula for n = 2, 4, 6. 

The £( 4 )-vertices with one external boson, (h) and (i), follow from 

J^ a [£ { %-mesons = ^ Tr(T a [^, <j>\) Tr( X ) (111) 
+ Tr T a (x d^<j) <f> + d^x<t> ~ <Pxd^ - 4> 8$ x) 
_^^Tr(T a [^,^]) 

and 

jiA^U-mesons = ^ Tr(T a [d^, <j>}) Tr(d„0 d v <j>) (112) 



+ 



^Tr(T a [^,0])Tr(^^) 

+ ^£ Tr T a {{d^ d v 4> - <p{d^, d v <t> d v 4>}) 

- ^{6Tr(T a [^,^]) Tr(x^) +TrT a ([^,<A 3 ] - 3 <f>[d^, <f>]<f>) Tr( X )} 

- ^{TrTaixld^,^ 2 + 2x<f> 2 d^<f>-2 0x{d^,<p 2 } + 2d^ X 4> 3 ) 

+ Tr T a (4 ^x4> + 3 [d^, 4>]x4> 2 + 3 (fx[d^, 4>] + 2 {d^, ^ 2 }x<P) 
+ Tr T a ( - 4 <f> 8^ M ~ 2 3 X - 2 d^> ^X + <t> 2 [d^>, 4>]x) } 

+ ^ ^{Tr T a { - 3 </>[3 M 0, + 2 c>^ - ^ <j> 2 8^) 

+ Tr T a (2 4> 2 [dfj,<j>, d u (j)} + 2 [d^, d v <\>\(\? - [d^ <p, d u (f> <j)) - [<j> d^, W])} 
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Finally, vertex (j) stems from and takes the following form in the notation of Q: 

J^ a [£ { %- meS ons = ^^T a [ x ,{d^,d^}} (113) 

iB 

+ -^pf Tr T a {2x4>df,4> x - 2x %4> <t> X + XX 9^ 4> 
- XX'pd^ + d^cpxx ~ 4>d^<Pxx) 



+ ^TY( X )TrT a [{^,0},x] 
iB 

+ ^ Ta ^ XX - <9 M XX<A + </>XX d^4> - XX 9^ 4>) 
iB 

+ ^2 Tr(x) Tr Ta^d^x - d^xfi + 4>xd^ ~ x9 lJ <t>4>) 
+ ^ {Tr(xx) Tt (Ta&, M) + 2Tr Mx, ^( T a 
+ ^TV(x9 M 0)Tr(r Q [0,x]) 
+ ^lT¥(x)Tr(x)Tr(T a [0,^]) 
+ ^ ItTaCft,^^^,^] - d^[8 v 
+ l J?*ff>df , TrT a ([d lt d( ) <f> t d v <t>] - [dudpfrdrf]) 

- l -j^d»TrT a (0^ X d v 4> - d v <f> X d^) 
-!|f^TV(T [^,^])TV(x) 
+ ^X d v TvTa ([d t ,cf>,{d u ct>,x}] ~ M,{W,X}] 



Note that in (0) we have given only the part of which is relevant for the K° form factor. 



We convinced ourselves by inserting the Feynman rule (112) into diagram (4) that only the terms 
containing B24 and B25 yield a nonzero contribution to the K° form factor. 
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